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Abstract. Spectra of the QSO pair Q0307-195A,B have 
been obtained in the Lya forest (3660-3930 A) and C iv 
(4720-4850 A) regions with a FWHM resolution between 
0.7 and 0.5 A. 46 hnes have been detected in the spectrum 
of object A while 36 in the spectrum of object B, of them 
29 and 20 were identified as Lya absorptions respectively. 
The present observations have been supplemented with 
data of comparable quality on other 7 QSO pairs available 
in the literature to give an enlarged sample of 217 Lya 
lines with rest equivalent width Wo > 0.3 A. The analysis 
of the hits (i.e. when an absorption line appears in both 
QSO spectra) and misses (i.e. when a line is seen in any of 
the QSO spectra, but no line is seen in the other) , carried 
out with an improved statistical approach, indicates that 
the absorbers have typically a large size|^: R — 362 
kpc, with 95% confidence hmits 298 < R < 426 kpc 
and R — 412 h^^ kpc, with 95% confidence limits 333 < 
R < 514 h^^ kpc for the radius of idealized spherical 
and disc geometries, respectively. The present data do not 
allow to establish any correlation of the typical inferred 
size with the proper separation or with the redshift of the 
pairs. 

The correlation between the observed equivalent 
widths of the absorbers in the adjacent lines of sight 
becomes poorer and poorer with increasing proper sep- 
aration. A disc geometry with a column density profile 
N(r) oc {r/Ro)~'^, 7 = 4, is found to reasonably repro- 
duce the data with Rq ~ 100 — 200 h^^ kpc, but also 
spherical clouds with the same column density profile and 
a power-law distribution of radii may give a satisfactory 
representation of the observations. 

Key words: intergalactic medium - quasars: absorption 
lines - quasars:individual(Q0307-195A,B) 
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* Based on observations collected at the European Southern 
Observatory, La Silla, Chile (ESO No. 149.B-0013). 

^ Throughout this paper we will assume: h = _ffo/(100 km 
s"^ Mpc"^), qo = 0.5. 



1. Introduction 

Closely separated QSO pairs and gravitationally lensed 
QSOs provide two or more adjacent lines of sight (LOS), 
which allow to sample the size and clustering of the ab- 
sorbers. 

The ray path separation of gravitationally lensed QSOs 
usually spans subgalactic scales; studies of these spectra 
have shown that the Lya clouds are much larger than a 
few kiloparsecs (Weymann & Foltz 1983; Foltz et al. 1984; 
Smette et al. 1992, 1995). Further results again pointing at 
very large sizes have been obtained for LOS separations 
of a few arcminutes or less (Shaver & Robertson 1983, 
SR83; Crotts 1989; Dinshaw et al. 1995; Dinshaw et al. 
1997, D97; Fang et al. 1996, FDCB; Crotts & Fang 1997, 
CF97) . The first indication in this sense has been obtained 
observing the QSO pair Q1343+2640A,B (with a proper 
sightline separation of 39—40 h^^ kpc at 2: ^ 2), indicating 
a cloud radius of 100 - 200 h"^ kpc (Bechtold et al. 1994; 
Dinshaw et al. 1994). QSO pairs with separations larger 
than a few arcminutes, corresponding to proper separa- 
tions up to megaparsecs, have the potential to study the 
large-scale spatial distribution of Lya or metal absorbers 
(Crotts 1985; Jakobsen et al. 1986; Tytler, Sandoval & 
Fan 1993; Elowitz, Green & Impey 1995; Dinshaw& Im- 
pey 1996; Williger et al. 1996). 

Multiple sightlines are also highly valuable in prob- 
ing the large scale structure. In particular, multiple, well- 
sampled sightlines can provide a stronger test for voids by 
placing more absorbers in a void-sized volume than could 
possibly be obtained along a single sightline. 

Furthermore, the effect of the UV radiation field of a 
foreground QSO on the absorption spectrum of a back- 
ground QSO provides an important way of testing the 
"proximity effect" interpretation (Bajtlik, Duncan & Os- 
triker 1988) of the "inverse effect" observed in single QSO 
sightlines. 

We present new observations of the QSO pair Q0307- 
195A,B discovered in an objective prism survey by 
MacAlpine and Feldman ( |l982| ) (A = UM 680 and B = 
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UM 681). Later on, the pair was observed by Shaver and 
Robertson ( |1983D who analyzed the absorption spectra of 
the two objects. Both QSOs are about of 19th magni- 
tude, their redshifts are Zom(A) = 2.1439 ± 0.0003 and 
0em(B) = 2.1217 ± 0.0003 (SR83), and they are separated 
by 56 arcseconds on the plane of the sky. The redshift dif- 
ference is small enough 2126 km s~^) that they could 
be members of the same cluster. 

The rest of the paper is organized as follows: in § 2, 
we present spectroscopic observations of Q0307-195A,B 
and describe the procedure of reduction and calibration 
of the spectra. In § 3 the metallic systems in the two QSO 
spectra are listed and briefly discussed. Section 4 describes 
the procedure to obtain an estimate for the radius of the 
Lya absorbers in the hypothesis of simple geometry. In 
subsection 4.1 we tackle the same issue on the basis of 
more realistic models for the absorbers. Finally, in § 5 we 
summarize our results. 

2. Observations and data reduction 

We obtained spectra of the QSO pair Q0307-195A,B with 
the blue arm of the EMMI spectrograph at the 3.5m New 
Technology Telescope of the European Southern Observa- 
tory in two runs in November 1994 and November 1995. 
The spectrograph was used with the holographic 3000 lines 
mm^^ grating (ESO # 11) and a 1024^, 24 /im pixel, 
Tektronix CCD (ESO # 31). Two or three spectra were 
obtained at each of three different central wavelengths: 
3690,3750 and 3850 A (Lya forest region) and one at 
4780 A (Civ region). Typical exposure times were 6000 
s/spectrum for a total exposure time of 56 900 s. In this 
instrument configuration one CCD pixel corresponds on 
average to 0.15 and 0.13 A in the UV and blue regions re- 
spectively. The typical slit width was 1.2 arcsec in width 
and 180 arcsec in length and it was aligned at P. A. 135°. 6 
to capture both QSOs. 

The spectra were extracted, wavelength calibrated, 
normalized and merged using standard procedures of the 
ESO MIDAS long slit reduction package. 

The final smoothed and merged spectra cover the 
ranges AA3660 - 3930 A and AA4720 - 4850 A. 

The resolution, as measured from the Thorium lines of 
the calibration spectra, extracted and treated in the same 
way as the QSOs spectra, is 0.7 and 0.5 A in the UV and 
blue range respectively. The s /n ratio of the final spectra 
varies from 25 at the peak of the Lya emissions of the two 
QSOs to 5 at the the shortest wavelengths. The relative 
accuracy of the wavelength scales in the spectra of the two 
QSOs is better than 20 km s~^. 

The normalized spectra are shown in Figs. |l|, ||. 

3. Detection and identification of lines 

The determination of the continuum in the QSO spec- 
trum is a critical step because it affects the measurement 



of the absorption line parameters. In the present work, we 
selected the portions of the spectrum free of strong ab- 
sorption hnes or artificial peaks (e.g. due to cosmic rays), 
i.e. where the RMS fluctuation about the mean becomes 
consistent with noise statistics, and then, we estimated 
the continuum level by spline-fitting these regions with 
quadratic polynomials. 

In order to detect the absorption lines we have searched 
for all the features deviating from the continuum more 
than 3.5 times the RMS noise. The lines have been fit- 
ted with Voigt profiles convolved with the instrumental 
spread function, making use of a minimization method of 
X^. This step has been performed within the ESO MIDAS 
lyman package (Fontana & Ballester 1995). The values 
of the wavelength, redshift z, and rest equivalent width 
Wo have been determined for isolated lines and individual 
components of blends. 

The number of components of each absorption feature 
is assumed to be the minimum required to give a reduced 

< 1 (corresponding to a confidence level P > 50%). 

The error on the observed equivalent widths was com- 
puted in such a way that it includes also the uncertainty 
in the determination of the level of the continuum and in 
the choice of the fitting interval for the lines. 

Due to the improved signal to noise ratio (s/n) and 
resolution, we observe more than twice the number of lines 
found by SR83 in the common wavelength interval. The 
lists of lines are presented in the Tables |l|, ^. 

3.1. Metal Systems 

We based the identification of metal lines on SR83. All the 
systems reported by them in the common wavelength in- 
terval are confirmed, in addition we found one new system 
in Q0307-195A and one possible new system in Q0307- 
195B (see 3.2 and 3.3). 

An identification programme, based on the method of 
Young et al. (1979), has been applied to the newly obser- 
ved lines that could not be attributed to systems already 
identified by SR83. Because of the limited wavelength cov- 
erage in the red region of the spectra, we were able to 
identify only one new Civ doublet. 

Four and five lines redward of the Lya remained 
unidentified in Q0307-195A and B respectively. 

The metal systems of each object are described in de- 
tail in the following section. 

3.2. Object Q 0307- 195 A 

The metal system at Zabs — 2.1229 - a 
Already observed by SR83. This absorption system 
presents a redshift difference of only 115 km with 
object B and of 77 km s~^ with system a of object 
B. SR83 argued that this common absorption is likely 
due to a very extended gaseous halo or disc, physically 
associated with Q0307-195B itself. 
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Metal 


In o. 








Identification 




I^^ystem 


1 


3661. 43± 0.04 


0.62 


1.88±0.04 


Si 111 (1206) 


2 


0348 


c 


2 


3663.31 ± 0.05 


0.16 


0.48±0.05 


H I (1216) 


2 


0134 




3 


3665. 80± 0.05 


0.17 


0.50±0.04 


H I (1216) 


2 


0155 




4 


3669. 36± 0.05 


0.33 


l.OOitO.lO 


H I (1216) 


2 


0184 




5 


3676.08± 0.06 


0.32 


0.97±0.06 


H I (1216) 


2 


0239 




6 


3689. 15± 0.02 


1.34 


4.07±0.03 


H T (1216) 


2 


0347 


c 


7 


3692. 82± 0.03 


0.19 


0.60±0.01 


H I (1216) 


2 


0377 




8 


3694. 16± 0.12 


0.23 


0.70±0.06 


H I (1216) 


2 


0388 




9 


3699. 64± 0.02 


0.26 


0.81±0.02 


H I (1216) 


2 


0433 




10 


3701.44± 0.04 


0.19 


0.59±0.03 


H I (1216) 


2 


0448 




11 


3704.01± 0.15 


0.20 


0.61±0.02 


H I (1216) 


2 


047 




12 


3704. 98± 0.08 


0.22 


0.68±0.02 


H I (1216) 


2 


0477 




13 


3716. 92± 0.08 


0.34 


1.03±0.04 


H T (1216) 


2 


0575 




14 


3719. 67± 0.03 


0.23 


0.71±0.03 


fl T (1216) 


2 


0598 




15 


3726.69± 0.04 


0.08 


0.26±0.01 


H T (1216) 


2 


0655 




16 


3730.29± 0.02 


0.41 


1.25±0.01 


H T (1216) 


2 


0685 




17 


3751. 39± 0.02 


0.54 


1.68±0.04 


H T (1216) 


2 


0859 




18 


3755. 94± 0.18 


0.08 


0.23±0.02 


H T (1216) 


2 


089 




19 


3757.60± 0.02 


0.42 


1.29±0.02 


H I (1216) 


2 


0909 




20 


3759.82± 0.03 


0.13 


0.40±0.03 


H I (1216) 


2 


0928 




21 


3766. 09± 0.13 


0.49 


1.52±0.02 


H I (1216) 


2 


097 




22 


3767. 22± 0.02 


0.43 


1.33±0.02 


H I (1216) 


2 


0989 


b 


23 


3767. 84± 0.06 


0.17 


0.54±0.02 


Si 111 (1206) 


2 


1229 


a 


24 


3770.53± 0.02 


0.15 


0.45±0.02 


H I (1216) 


2 


1016 




25 


3774.58± 0.04 


0.05 


0.16±0.01 


H T (1216) 


2 


1049 




26 


3776. 67± 0.02 


0.16 


0.50±0.03 


H I (1216) 


2 


1066 




27 


3779. 91± 0.02 


0.13 


0.41±0.03 


H I (1216) 


2 


1093 




28 


3796. 55± 0.01 


0.67 


2.10±0.05 


TI T (1216) 


2 


1230 


a 


29 


3800.98± 0.01 


0.43 


1.33±0.03 


H T (1216) 


2 


1266 




30 


3805. 15± 0.07 


0.03 


O.lliO.Ol 


H I (1216) 


2 


1301 




31 


3817. 25± 0.01 


0.12 


0.37±0.02 


H I (1216) 


2 


1400 




32 


3825.47± 0.01 


0.17 


0.50±0.02 


Si 11 (1260) 


2 


0351 


c 


33 


3838. 96± 0.02 


0.06 


0.19±0.02 


N V (1238) 


2 


09888 


b 


34 


3851.31± 0.02 


0.04 


0.13±0.02 


N V (1242) 


2 


09888 


b 


35 


3867. 99± 0.07 




0.14±0.02 


not ident. 








36 


3868. 81± 0.01 


0.14 


0.43±0.01 


N V (1238) 


2 


12298 


a 


37 


3881. 25± 0.01 


0.13 


0.39±0.02 


N V (1242) 


2 


12298 


a 


38 


3897.72± 0.06 




0.51±0.02 


not ident. 








39 


3911. 83± 0.04 


0.36 


0.92±0.05 


C TV (1548) 


1 


5267 


d 


40 


3916. 17± 0.04 




0.23±0.02 


not ident. 








41 


3918. 34± 0.04 


0.30 


0.77±0.05 


C TV (1550) 


1 


5267 


d 


42 


3927.40± 0.02 




0.58±0.04 


not ident. 








43 


4797. 22± 0.02 


0.24 


0.74±0.03 


C TV (1548) 


2 


0986 


b 


44 


4805.20± 0.02 


0.19 


0.59±0.02 


C TV (1550) 


2 


0986 


b 


45 


4834. 71± 0.02 


0.33 


1.03±0.08 


C TV (1548) 


2 


1228 


a 


46 


4842. 75± 0.02 


0.23 


0.72±0.04 


C TV (1550) 


2 


1228 


a 
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Table 2. Absorption lines in the spectrum of Q0307-195B 



A (Obs.) Wa (rest) Wa (obs) 



1\ (J. 




/A) 










i A ) 




IdcDtificQitiOD 


2 


lbs 


SystGni 


1 


3660 


76± 





47 


0.72 


2 


18±0 


18 


Si TJJ (1206) 


2 


034 


b 


2 


3662 


71± 





15 


0.40 


1 


20±0 


17 


H I (1216) 


2 


0129 




3 


3669 


45± 





10 


0.54 


1 


64 ±0 


04 


H I (1216) 


2 


0185 




4 


3672 


85± 





11 


0.32 





96±0 


07 


H I (1216) 


2 


0213 




5 


3677 


81± 





08 


0.27 





83±0 


03 


H I (1216) 


2 


0253 




6 


3684 


11± 





08 


1.16 


3 


52±0 


04 


H I (1216) 


2 


0305 




7 


3687 


80± 





11 


1.19 


3 


62±0 


04 


H I (1216) 


2 


0335 


b 


8 


3693 


10± 





15 


0.54 


1 


65±0 


13 


TT 1 (1216) 


2 


0379 




9 


3696 


62± 





05 


0.14 





'l'l±0 


02 


IT I (1216) 


2 


0'I08 




10 


3700 


04± 





04 


1.25 


3 


80±0 


04 


H I (1216) 


2 


0436 




11 


3703 


81± 





06 


0.28 





85±0 


03 


H 1 (1216) 


2 


0467 




12 


3721 


32± 





04 


0.99 


3 


02±0 


09 


H 1 (1216) 


2 


0611 




13 


3723 


77± 





04 


0.49 


1 


51±0 


09 


H I (1216) 


2 


0631 




14 


3730 


48± 





05 


0.17 





53±0 


03 


H I (1216) 


2 


0687 




15 


3735 


50± 





05 


0.23 





72±0 


07 


H 1 (1216) 


2 


0728 




16 


3767 


14± 





01 


0.64 


2 


00±0 


04 


Si III (1206) 


2 


1224 


a 


17 


3769 


23± 





10 


0.08 





25±0 


04 


H 1 (1216) 


2 


1005 




18 


3770 


50± 





15 


0.10 





32±0 


04 


H 1 (1216) 


2 


1016 




19 


3775 


53± 





04 


0.04 





14±0 


03 


H I (1216) 


2 


1057 




20 


3780 


09± 





01 


0.16 





49±0 


02 


H I (1216) 


2 


1095 




21 


3782 


72± 





02 


0.53 


1 


66±0 


05 


H I (1216) 


2 


1116 




22 


3795 


52± 





01 


0.60 


1 


86±0 


05 


II I (1216) 


2 


1222 


a 


23 


3821 


93± 





01 


0.19 





57±0 


03 


Si II (1260) 


2 


0323 


b 


24 


3847 


53± 





22 







25±0 


04 


not ident. 








25 


3880 


10± 





18 







19±0 


02 


not ident. 








26 


3885 


27± 





06 


0.08 





22±0 


05 


Si IV (1393) 


1 


7876 


c 


27 


3886 


88± 





0-1 


0.81 


2 


27±0 


05 


Si IV (1393) 


1 


7888 


c 


28 


3904 


40± 





15 







22±0 


04 


not ident. 








29 


3909 


60± 





05 


0.15 





38±0 


05 


C IV (1548) 


1 


5253 


d 


30 


3910 


40± 





06 


0.05 





14±0 


05 


Si IV (1402) 


1 


7876 


c 


31 


3912 


02± 





04 


0.57 


1 


60±0 


08 


Si IV (1402) 


1 


7888 


c 


32 


3916 


10± 





05 


0.08 





21±0 


06 


C IV (1550) 


1 


5253 


d 


33 


3922 


22± 





04 







43±0 


02 


not ident. 








34 


3925 


67± 





04 







40±0 


07 


not ident. 








35 


4833 


32± 





02 


0.51 


1 


59±0 


11 


C IV (1548) 


2 


1219 


a 


36 


4841 


36± 





02 


0.44 


1 


37±0 


05 


C IV (1550) 


2 


1219 


a 
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Fig. 1. Spectra of Q0307-195A,B obtained with the NTT EMMI instrument as a function of vacuum hehocentric 
wavelength. The dashed line shows the la error in the flux. Tick marks indicate absorption lines detected at or above 
the 3.5a confidence level, upper and lower tick marks are referred to Lya and metal absorption lines respectively 



The metal system at Zabs = 2.0988 - b 
This is a new system. We identified the C iv and N v 
doublets in the red portion of the spectrum, and the 
Lya line blended with the Sim 1207 of the previous 
system. 

The metal systems at z^hs = 2.03497 - c and ^abs = 
1.5267 - d 

Already observed by SR83. No new line has been 
added. 

3.3. Object Q0S07-195B 

The metal systems at z^hs = 2.1221 - a and Zabs = 2.0333 
- b 

Already observed by SR83. No new line has been 
added. 

The metal system at ^;abs = 1-7882 - c 
Already observed by SR83. The identification is based 
only on the Si IV doublet, in the Lya forest a possible 
Cii can be identified but it is blended with a strong 
Lya line. The Si IV doublet shows velocity structure 
and it was fitted with two components. 

The metal system at Za,bs = 1-5253 - d 
The identification of this possible new system is based 
only on the CiV doublet found redward of the Lya 



emission line. This system differs in redshift by only 
166 km s~^ from system d in object A. 

All the lines in the Lya forest not identified as metals 
are assumed to be Lya. 

4. The size of the Lya absorbers 

The main aim of this work is to determine the transverse 
size of the Lya absorbers using the Lya forests of the 
two adjacent LOSs of the pair. In order to avoid spurious 
effects due to the "proximity effect", only the Lya lines 
with velocity separation greater than 5000 km s~^ from 
the lower QSO redshift in the pair were considered in the 
sample. 

The spectra of the observed pairs show a number of 
hits, J\f^, and of misses, Af^. A hit occurs when an ab- 
sorption line above a given threshold in rest equivalent 
width. Wo, appears in both QSO spectra, with a velocity 
difference less than a limit, Av. A miss occurs when a line 
is seen in any of the QSO spectra, but no line above the 
Wo threshold is seen in the other. 

A value Av = 200 km s~-^ has been adopted as the 
most suitable velocity window, being significantly larger 
than the resolution (~ 40 km s^^) and close to the clus- 
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-d 

N 




0.5 



4800 



4850 




4800 
Wavelength (A) 

Fig. 2. Same as Fig. ^ for the C iv region 

Table 3. Lya lines that give rise to hits in the spectrum of Q0307-195A,B. 



4850 



Aa (Obs.) Wa (rest) Ab (Obs.) 



No. 


(A) 


(A) 


(A) 


(A) 


(km s ^) 


1* 


3669.36 


0.33 


3669.45 


0.54 


-7±9 


2 


3676.08 


0.32 


3677.81 


0.27 


-141±8 


3* 


3689.15 


1.34 


3687.80 


1.19 


110±9 


4 


3694.16 


0.23 


3693.10 


0.54 


86±16 


5 


3699.64 


0.26 


3700.04 


1.25 


-32±4 


6 


3704.01 


0.20 


3703.81 


0.28 


16±13 


6 


3704.98 


0.22 


3703.81 


0.28 


95±8 


* Hit 


in the case Wo 


> 0.3 









tering scale observed for Lya absorbers (Cristiani et al. 
1997). As a consequence, all the Lya lines in the sam- 
ple whose separation was less than 200 km s^^ have been 
merged into a single line with wavelength given by the 
average of the wavelengths and equivalent width equal to 
the sum of the equivalent widths of the "parent" lines. 

In Table ^, the hits observed in the spectra of Q0307- 
195A,B are reported together with the velocity differences 
between the lines. 

To constrain cloud sizes, it can be assumed that a 
binomial random process produces the observed number 
of hits and misses; if the probability for a hit is given 
by VE*, the likelihood function of this binomial process is 
/;(A/h°,A/'°|^') = ^-^h (1 - ^')^» (FDCB). 



5* is the probability that both ray paths intersect the 
cloud, given that at least one ray path does. ^E" is related to 
the probability <i>, that one given LOS intersects the cloud 
given that the other adjacent ray path already does, by 
the formula: 



* = $/(2-$). 



(1) 



The functional form of the probability $ depends 
on the assumptions for the geometrical shape of the 
absorbers. In the hypothesis of single radius, spherical 
clouds, with the definition X = S/2R, where S is the 
proper separation of ray paths and R is the cloud radius, 
it is found (McGiU 1990): 



$ = (2/7r) [arccosX -X{1- X 



2U/21 



for X <l, 



(2) 
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Table 4. Parameters of the QSO pairs forming our enlarged sample and Lya cloud radius estimates in the hypothesis 
of spherical and disc geometry 



Spheres Discs 





Angular 


Lya z 


S, Proper 


Wn 






95% Confid. 


Modal 


95% Confid. 


Modal 


QSO Pair 


Separ. 


Range 


Separation 


ciil off 


K 


M" 

* m 


Inlerval in /? 


Radius 


Inlerval in /? 


Radins 




(arcsec) 




(/(-I kpc) 


(A) 






(h-^ kpc) 


kpc) 


(h-^ kpc) 


(h-^ kpc) 


1343+2640A/B1 
1026-0045A/B^ 


9.5 


1.756-1.979 


39-40 


0.4 


6 


1 


50-2705 


80 


70-3000 


117 


36 


0.833-1.398 


149-155 


0.2 


5 


15 


98-268 


137 


117-385 


179 










0.3 


4 


12 


98-305 


139 


118-443 


183 


0307-195A/B3 


56 


2.011-2.070 


228-229 


0.2 
0.3 


6 
2 


7 

9 


10-294 
< 269 


146 


11-406 
< 365 


171 


0107-025A/B'* 


86 


0.481-0.892 


301-361 


0.3 


5 


6 


291-1776 


475 


379-2636 


663 


1517+2356/57= 
1623+2651A/B= 


102 


1.390-1.787 


427-438 


0.4 


3 


21 


34-460 


277 


39-618 


330 


127 


2.025-2.409 


497-519 


0.4 


5 


20 


99-604 


348 


113-828 


425 


1623+2651A/53= 
1623+2653/51B= 


147 


1.971-2.409 


575-604 


0.4 


7 


18 


356-883 


481 


414-1247 


619 


177 


2.025-2.468 


687-723 


0.4 


5 


32 


27-600 


418 


30-773 


482 



1 Crotts et al. 1994 

2 Pctitjean et al. 1998 
^ This paper 

* Dinshaw et al. 1997 
= Crotts & Fang 1997 



and i> = otherwise. 

As a second model, clouds can be idealized as circu- 
lar discs, with a given radius R, much larger then their 
thickness, and an observed inclination angle 6. The prob- 
ability that one ray path intersects the disc within the 
angle 9 ^ 9 + d9, given that the other ray path already 
does, is (McGill 1990): 



$(0) = {cos9/tt)1 arccos 



1- 



X 
cos 9 



X 
cos 9 

1/2 



X 
COS 9 



for X < COS 9, 



(3) 



and <I> = otherwise. By integrating over 9 for randomly 
oriented discs, the probability is then given by: 



r/2 

/ <^{9)d9. 

-Tr/2 



(4) 



In practice, from the observations we determine the 

numbers of hits and misses, J\f^ and J\f^, and then wc 
correct them for the accidental hits due to random velocity 
matches. 

Two different procedures have been devised to com- 
pute the value of Af^: in the first, the velocity separation 
between all the line pairs is computed and, when a given 
line combines to give more than one hit, only the one with 
the smallest velocity separation is taken into account, and 
the others are neglected. If, for example, lines 1 and 2 of 
spectrum A combine with lines 1 and 2 of spectrum B 
forming the couples 1-1, 1-2 and 2-2, and the couple 1-2 



has the smallest velocity separation, we count just one hit. 
The second method takes the minimum value between the 
"A vs. B" hits and the "B vs. A" hits. 

The two methods turned out to yield the same result 
in all the investigated cases. 

The number of observed misses is given by the sum of 
the two Af^'s, "A vs. B" and "B vs. A". 

To take into account the possibility of accidental hits 
due to random velocity matches, we have used simulations 
of Poissonian distributions of absorption lines in the two 
QSO LOSs, with a number of lines in each spectrum equal 
to the observed one. 

For each simulation we have computed: 

- the number of accidental hits, A/'r, applying the same 

procedures used for the real spectra; 

- the corrected number of hits, Mh, (misses, Afm) by sub- 

tracting A/'r from (adding 2 A/'r to) the observed value. 

Eventually, the probability distribution of the number 
of corrected hits C(A/'h) was evaluated. The probability 
C(0) takes into account all the occurrences with A/'r > 

For each pair of values of corrected hits and misses, 
given the expressions for the probability distributions $ 
and 4*, it is possible to apply Bayes's Theorem (Press 1989, 
FDCB) to yield the a posteriori probability density for the 
radius R, 



V{R,Afh,Afm) = 



r(A4,A4r|^)/(J?) 

/o°°£(A4,A4r|*)/(i?)rfi?' 



(5) 



At variance with FDCB, we have chosen as a prior 

distribution, f{R), a uniform distribution for R between 
and i?max) and not a uniform distribution of ^'(i?)- 
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Fig. 3. Probability distribution V{R) as given by eq. (8) 
for spherical absorbers (solid curve) and randomly in- 
clined discs [dotted curve) as functions of cloud radius for 
all the QSO pairs considered in the paper. The plateau 
at separation below S/2 arises because there is a non- 
vanishing probability that the number of corrected hits 
be (C(0) 7^ 0) 



It is worth stressing that a uniform prior distribution 

for the probability 4* introduces a cutoff at high and low 
R that forces the 95 % confidence intervals to be narrow 
also in situations with a very low number of hits or misses. 
In particular, a spurious dependence of R on the proper 
separation S is induced in the case of very few hits, which 
is the common situation at large separations, because the 
probability '^{R) actually depends on S. Since there is 
no particular justification to assume a uniform \E'(-R), we 
have chosen a uniform distribution for R which does not 
introduce biases in the results. 
Equation (5) then becomes 



V{R,M^,M^ 



/:(AAh,AA^|*) 



(6) 



We have adopted as i?max a value of 3 Mpc which 

represent the mean free path between two absorbers along 
the line of sight, as obtained from the line distribution 
d'^n/dzdNm (Giallongo et al. 1996) at a redshift ^ ~ 2, 
in the column density range 14 < log A''hi < 17. 

In the case of discs, the maximum value is increased 
by a factor ^/2, yielding a value i?max — 4/i~^ Mpc. 

It has to be noted that the value of this maximum 
radius does not influence in practice the resulting modal 



Fig. 4. Inferred cloud radius R (and 95 % confidence in- 
terval) from a model assuming single-radius, unclustered 
spherical clouds, as a function of QSO pair sightline sepa- 
ration S. The dashed line shows the most probable value 
for R and the dotted lines the 95 % confidence interval 
(see text) 

radii and confidence intervals, because the probability dis- 
tributions for R usually go to zero well before this upper 
limit. i?niax plays a significant role only if the number of 
misses is extremely low (1 or 0). 

In the case of no hits, recipe (6) gives an incorrect 
result, in fact, the probability $ becomes zero for R < S/2 
by definition. In this extreme occurrence we have then to 
write: 



V{R,Mi„M^) = Ax 



1 if i? < S/2 

C{Afh,Afm\^) if i?> 5/2 



(7) 



where, A is the proper normalization factor which makes 

/(f'"^''7'(i?,Mi,AA^) = 1. 

The global V{R) for a given QSO pair is finally ob- 
tained by the sum: 



V{R) = ^V{R,Mh,Mm) X C(AAh). 



(8) 



In order to investigate the possible correlations 
of the radius R with the LOS separation and with 
the redshift, we have applied our analysis to 7 
other data samples of comparable quality, pub- 
hshed in the literature: Q1343+2640A,B (Crotts et 
al. 1994), Q1026-0045A,B (Petitjean et al. 1998), 
Q0107-025A,B (D97), Q1517-h2356,1517-h2357 (CF97), 
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Q0107-025A/B 
<z>=0,7 



Q1343+2640A/B 
<z>=l,9 



Q 1633 +265 lA/B 
<z>=3.2 



Q0307-195A/B 
<z>=2.1 



Q1517-2357/56 
<z> = l,6 



Q10a6-0045A/B 



Q1623-I-2651A/53 
<z> = g.2 



Qi623+2653/51E 
<z>=2.3 



200 400 600 

Separation S (h ' kpc) 



Fig. 5. Same as Fig. ^ but for single-radius, unclustered 
thin discs 



Q1623+2651A,1623+2653,1623+2651B (CF97). The 
QSO pair Q 1343+2640 A, B has been observed also by 
Dinshaw et al. ( 1994 ), we have chosen to use Crotts et al. 



( 1994 ) line list because of the larger wavelength coverage 
and larger number of identified lines. 

The relevant information are displayed in Table ^. 
Columns 6 and 7 report the observed number of hits and 
misses, respectively. 

The most probable number of accidental hits for each 
QSO pair is: for 1343+2640A,B; for 1026-0045A,B 
(for both Wo = 0.2 and Wo = 0.3); 2 for 0307-195A,B 
(Wo = 0.3); 6 for 0307-195A,B {Wo = 0.2); for 0107- 
025A,B; 1 for 15174-2356/57; 2 for 1623+2651A,B; 2 for 
1623+2651A/53 and 5 for 1623+2653,51B. 

In the actual computation of the cloud sizes the aver- 
age value of the proper separation has been adopted. 

In Fig. ^ the probability distributions given by equa- 
tion (8) for all the QSO pairs presented in the paper are 
plotted for the two cases of spherical and disc geometry. 

In Fig. ^ and Fig H the values estimated for the modal 
radii are plotted as a function of S. We consider that if 
the typical radius R^o < S/2, where R^q is defined as 



V{R) = 0.5, 



(9) 



only an upper limit can be set to the dimension of the ab- 
sorber. This is the case for Q0307-195A,B (with Wo > 0.3 
A) for which R^q ~ 90 h^^ kpc < 114 h^^ kpc for the 
spherical model and R^q ~ 107 h^^ kpc < 114 kpc 
for the disc model. The dashed line represents the most 



probable modal radius, R — 362 kpc, with 95% con- 
fidence limits 298 < R < 426 /i^^ kpc (spherical ab- 
sorbers) and R — 412 h^^ kpc, with 95% confidence limits 
333 < i? < 514 kpc (disc). These values correspond 
to the peak of the total probability density obtained as 
the product of the probability densities of the considered 
QSO pairs. 

A plot of the modal radius as a function of the average 
redshift of the Lya forest of the corresponding QSO pair 
does not show any correlation between the two, at variance 
with the result obtained by Dinshaw et al. ( 1998| ). 

The quality of the present data on the pair Q0307- 
195A,B and of the data by Petitjean et al. (|1998D for the 



pair Q1026-0045A,B allows defining a complete sample of 
lines down to an equivalent width Wo — 0.2 A. The result- 
ing number of hits and misses, the radius and confidence 
interval are reported in Table ^. The corresponding values 
of the global modal radius and 95 % confidence interval are 
R = 350 kpc and 276 < i? < 417 h'^ kpc (sphere), 
R = 408 h-^ kpc and 325 < i? < 501 /i"! kpc (disc). 

Until fairly recently it has been conjectured that the 
typical transverse sizes of the Lya were of the order of a 
few tens of kpc. The present analysis indicates about one 
order of magnitude larger dimensions at z ~ 2, implying 
correspondingly larger ionizations and masses. With the 
presently inferred sizes and the typical UV background 
conditions at z ~ 2 (Giallongo et al. 1996), optically thin 
absorbers of log(iVHi) '-^14—15 are expected to show a 
remarkably small fraction of neutral Hydrogen, roughly 
2Q-5.5 _ 10^^. In this way the mean intergalactic density 
contributed by these clouds (Hu et al. 1995), especially in 
the spherical case, is close to conflict with the baryon limit 
from nucleosynthesis Of, < 0.015/i~^ (Walker et al. 1991). 
This potential incongruity can however be alleviated by 
assuming a disc or other geometries. 



.1. Towards 



realistic models 



The approach described in the previous section, although 
customarily used to determine the size of the Lya ab- 
sorbers, is expected to give just a rough idea of the dimen- 
sion of these structures, due to the simplistic hypotheses 
made on their geometry. 

In particular, the information on the equivalent width 
of the lines are not taken into account, since no density 
distribution is assumed inside the absorbers. Besides, all 
the clouds are assumed to have the same radius. 

In the literature (Smette et al. 1992; Smette et al. 1995; 
Charlton et al. 1997; Dinshaw et al. 1997) this problem has 
been tackled adopting a statistical technique that utilizes 
the information about the equivalent width of the lines. 

Dinshaw et al. (1997) adopted a power-law column 
density profile for the absorbers, N{r) = N\i^{r / Rq)^'' , 
where Mim = 1-26 x lO^"' cm^^ for 6 = 35 km s^^ is the 
limiting column density of the sample, Rq is the radius of 
the absorber at = A^um and 7 = 4, and three simple 
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Fig. 6. Equivalent widths of the hits and misses for the 8 QSO pairs in our sample plus the data by Smette et al. (1995), 
they are grouped as follows: (a) HE1104-1805; (b) Q1343+2640A,B; (c) Q1026-0045A,B Q0307-195A,B and Q0107- 
025A,B; (d) Q1517+2356, 1517+2357 Q1623+2651A,B and Q1623+265lA,1623+2653; (e) Q1623+2653,1623+2651B. 
The equivalent widths are arranged such that Wi = meux{WA, Wb) and W2 = min(iyyi, Wb)- The horizontal dashed 
lines represent the rest equivalent width threshold: Wo > 0.2 for (a), Wo > 0.3 for (b), (c) and Wo > 0.4 for (d), (e). 
The error bars in the upper right corners of the plots represents the average 1 a errors for the equivalent width 



geometries with fixed radius: spheres, randomly inclined 
discs, and pseudo-filamentary structures (approximated as 
discs with fixed inclination cosi = 0.2). 

The column density distribution along paired lines of 
sight for each of the geometric models considered is simu- 
lated by means of the Monte Carlo techniques outlined 
by Smette et al. (1992). This distribution can be con- 
verted into the corresponding equivalent width distribu- 
tion adopting an approximation to the curve of growth 
(Chernomordik & Ozernoy 1993). 

In Fig. 1^ the equivalent widths of the hits and misses 
are plotted for an enlarged sample of QSO pairs formed by 
our previously described 8 cases plus the lensed QSO HE 



1104-1805A,B observed by Smette et al. ( |1995| ). The pairs 
are grouped according to their typical proper separation. 
It should be noted that it is now certain that the dou- 
ble QSO HE 1104-1805A,B is actually a gravitational lens 
system, as the lensing galaxy has been detected (Courbin 
et al. 1998; Remy et al. 1998). The exact values of the sep- 
aration between the two light beams depend on the still 
badly determined lens redshift. The values we reported in 
Fig. ^ have been estimated using Fig. 6 of the paper by 



Smette et al. ( |1995| ), for the redshift zicns ~ 1.32, favoured 
by Remy et al. (1998). 

It is clear from the figure that the fraction of hits de- 
creases as the proper separation increases, while, at the 
same time, the correlation between the equivalent widths 
of the absorbers in the two LOS becomes poorer and 
poorer. 

The present observations can be compared with the 
simple model predictions of D97 (their Fig. 12). In particu- 
lar, our Fig. ||(a) shows reasonable agreement with the pre- 
diction of a disc model with i?o ~ 7 5* that is Rq ^ 90 
kpc; panel (b) agrees with the case Rq — 2 — 3 S* for 
both filament and disc model, that is i?o — 80 — 120 
kpc; panel (c) can be explained by filaments and disc with 
Rq — S which implies Rq ~ 220 kpc and, finally, 
panels (d) and (e) agrees with filament and disc models 
with i?o < "S* or i?o < 500 - 700 h'^ kpc. Altogether a 
disc geometry with Rq ^ 100 — 200 h^^ kpc seems to be 
favored. 

It should be noted, however, that the number of de- 
grees of freedom is not the same in the three models: once 
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a given radius is chosen in the disc idealization the two 
orientation angles can vary at random, for the filament 
geometry one angle is fixed, and for the sphere no angu- 
lar variable exists to define an orientation. In this way, 
it is obviously easier for the disc model to reproduce the 
real data that show a substantial amount of scatter due 
to measurement errors and possibly to a non-unique value 
of the cloud size. 

In fact, by introducing one more degree of freedom 
in the spherical cloud model: i.e. allowing the radius R 
to vary according to one of the two distribution: (l)expo- 
nential, and (2)power law with indices a = —4, —3, 5, 6, 
it is found that a power law distribution with a = 3 pro- 
vides an equally reasonable representation of the observed 
equivalent width distribution (D97). 

5. Conclusions 

We have presented new spectra of the quasar pair Q0307- 
195A,B. 

The number of detected lines is almost doubled with 
respect to the previous study of this pair (Shaver & 
Robertson 1983). One new metal system has been identi- 
fied in object A on the basis of the C iv doublet and one 
possible system is found also in object B. 

In order to study the shape and dimension of the 
Lya absorbers, we have added to the observations of the 
present QSO pair, data of comparable resolution taken 
from the literature to create a sample of 8 QSO pairs. We 
looked for hits and misses in the pair spectra and com- 
puted the probability distribution for a given radius R in 
the case of spherical and disc geometry basing our sta- 
tistical analysis on the approach outlined by Fang et al. 
( 119961 ). We modified some of their assumptions in order to 
deal correctly with the extreme cases of very low numbers 
(1 or 0) of hits or misses and avoid any spurious depen- 
dence on the separation between lines of sight. 

Our results can be summarized as follows: 

- the total probability density obtained as the product of 

the probability densities of the considered QSO pairs, 
gives a modal radius R — 362 h^^ kpc, with 95% con- 
fidence limits 298 < R < 426 h^^ kpc (spherical ab- 
sorbers) and R — 412 h^^ kpc, with 95% confidence 
limits 333 < i? < 514 h'^ kpc (disc). 

- No significant correlation is detected between the modal 

radius and the mean proper separation of the pair, nor 
between the modal radius and the mean redshift of the 
Lya forest region. 

- The comparison of the distribution of the equivalent 

widths of the observed coincident lines with the pre- 
dictions of Monte Carlo simulations for three simple 
geometrical models, spheres, discs and filaments with 
a power-law column density profile N{r) cx {r/Ro)^'' 
(Dinshaw et al. 1997), seems to indicate a disc shape 
with Rq ^ 100 — 200 kpc and 7 '--^ 4 as the most 



suitable. This could however be simply the result of the 
different number of parameters in the three models. If, 
for example, a further degree of freedom is allowed to 
the model based on spheres in the form of a power-law 
distribution of radii, a satisfactory consistency with 
the data is obtained. 
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